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Abstract 

We propose a scheme to investigate the interplay between Kondo-exchange interaction and quan¬ 
tum spin Hall effect with ultracold fermionic alkaline-earth atoms trapped in two-dimensional op¬ 
tical lattices using ultracold collision and laser-assisted tunneling. In the strong Kondo-coupling 
regime, though the loop trajectory of the mobile atom disappears, collective dynamics of an atom 
pair in two clock states can exhibit an unexpected spin-dependent cyclotron orbit in a plaquette, 
realizing the quantum spin Hall effect of the Kondo singlet. We demonstrate that the collective 
cyclotron dynamics of the spin-zero Kondo singlet is governed by an effective Harper-Hofstadter 
model in addition to second-order diagonal tunneling. 
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I. INTRODUCTION 


In the milestone work by Jun Kondo, the resistance minimum at nonzero temperature 
in a metal was explained by the scattering of conduction electrons by localized magnetic 


impurities 


Q. 


Since that, Kondo effect has been considered as a primary mechanism in 


heavy-fermion systems Q|. After the description of this effect in the ground state of the 
Kondo lattice model (KLM) was formulated, the formation of Kondo singlet between mobile 
electron and localized spin is found to perform a key function [3 6]. Thanks to the rapid 
development of laser cooling and trapping techniques, neutral atoms have been considered 
as an ideal platform for simulating complicated phenomena in condensed matter physics Q. 
For example, quantum mixtures of alkali atoms are expected to display the characteristics of 
Kondo-correlated state Q and multi-orbital effect foj . In particular, fermionic alkaline-earth 
atoms (AEAs) in optical lattices possess unique properties for operating optical atomic clocks 
with unprecedented precision 10, [111] and provide novel insights into the physics of strongly 


correlated transition-metal oxides and heavy-fermion materials as quantum simulators 


Q 


12 


16|, where the interorbita. 


plays an essential role 17 


spin-exchange interaction between 1 S 0 and 3 P 0 clock-state atoms 


19]. 


Active interest has recently been focused on operating neutral atoms as charged parti¬ 
cles by engineering artificial gauge potentials 


20 


24], As a significant breakthrough, spin 


Hall effect has been observed in a Bose-Einstein condensate, leading to the realization of 
a cold-atom spin transistor 25], where the internal and external states of atoms are cou¬ 
pled by a two-photon Raman process. For trapped atoms in optical lattices, tunable gauge 


potentials have been implemented by periodic driving 


26 ] and “Zeeman lattice” [27] tech¬ 


niques. Furthermore, the Harper-Hofstadter Hamiltonian 


28 


realized with ultracold atoms in optical lattices 


30 


29| has been experimentally 


31] using laser-assisted tunneling (LAT) 


321 ] . leading 


effect 


33 


o a spin-dependent magnetic field and realizing the quantum spin Hall (QSH) 
34], Moreover, since the interplay between interaction and QSH effect leads to 

sive bosons in optical lattices [35] 

, 4] in graphene or quantum wells 


compelling novel physics, both cyclotron motion of repu 
and QSH insulators with Kondo-exchange interaction [3 


368-42] have been extensively discussed. 


Motivated by these works, we present a proposal to study the interplay between Kondo- 
exchange interaction and QSH effect with ultracold fermionic AEAs trapped in two- 
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dimensional optical lattices using ultracold collision and LAT. The underlying physics 
is probed by the cyclotron dynamics of a mobile atom and a Kondo singlet in a plaquette. 
We find that the cyclotron orbit of mobile atom is damped with nonzero Kondo coupling. 
More interestingly, at strong Kondo coupling, we predict a striking and unexpected phe¬ 
nomenon that the spin-zero Kondo singlet can behave as a spin-half particle exposed to a 
perpendicular spin-dependent magnetic flux. We demonstrate that the spin state of local¬ 
ized background atoms and second-order diagonal tunneling have nontrivial effects on the 
spin-dependent cyclotron dynamics of the Kondo singlet. 


II. MODEL 


We consider fermionic AEAs in 1 S 0 (g) and 3 P 0 (e) clock states independently trapped 
in two-dimensional optical lattices of the same periodicity 43] as sketched in Fig. [[] where 
/ -I) denotes nuclear spin states 1 1, ±mj) [12 14]. The \\ / \) g atom (conduction electron) 
receives opposite frequency detuning 5cu = ±(cu 2 — oq) = ±A and momentum transfer 5 k = 
±(&2 — k \) from two far-detuned Raman beams and tunnels spin-dependently in the optical 
lattice tilted by a magnetic held gradient A along x 33]. In addition, the Mott insulator 
background of e atoms (localized spins) interacts with the g atom via onsite interorbital spin- 
exchange interaction V ex oc ( a~ g — a+) f d 3 rw g (r)w^(r), where af g is the scattering lengths 
for two atoms in |±) = -^=(| ge) ± | eg)) and w g ( e )(r) denotes the wavefunction for the g (e) 


atom 


HE, 


3l|. By averaging out rapidly oscillating terms in a rotating frame 3(1,13.1], we 


obtain the time-independent Harper-Kondo Hamiltonian (see Appendix A): 


H = H t + H k 

= ~ ^ 1 {Ke < ^ rn ’ n ' 0 ‘c\ n+1 ^ nol C m ,n,a + '^ C m,n+l,a C m,n,a + h.C.) 

m,n,a 

+ Vex y ] s m,n ’ ® m,ni (^) 


where s r m n = \ C L,n,a^af}C m , n ,g and S( n n = fL,n,a^ocgfm,n,p denote the spin operators 

a, (3 a, (3 

of the g and e atoms at (m,n) (a, (3 G {t, 1} and a is the vector of Pauli matrices). c\ nna 
( fLna ) creates a g (e) atom in the spin state |ct) at (m,n). is the spin-dependent 

complex tunneling amplitude of the mobile atom in the x direction induced by LAT, while 
J is the real tunneling amplitude of the mobile atom in the y direction. The first term 
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FIG. 1. (Color online) (a) Laser-assisted tunneling and Kondo coupling for |f) (top) and |f,) 
(bottom) g atoms. Though normal tunneling in the x direction is suppressed by a magnetic field 
gradient A for two g atom spin states, two Raman beams with frequency detuning 5u = U 2 —UJ 1 = A 
and momentum transfer 5k = k -2 — k\ restore the resonant tunneling with a spin-dependent complex 
amplitude Ke ±l ^ >rn ’ ri . (b) Schematic picture of effective magnetic flux in a plaquette with four sites 
A, B, C , and D , where the cyclotron orbits of the g atom that realize the spin-dependent effective 
magnetic flux ±<f> are coupled to a plaquette of localized e spins via onsite interorbital spin-exchange 
interaction (Kondo coupling) V ex . Ht denotes the Harper Hamiltonian. 


in Eq. (JT|) , i.e. the Harper Hamiltonian Ht 28|, 29], describes the g atom in a spin- 
dependent magnetic field arising from spatially varying phase 0 m , n ,t/! — ±0 m ,n = ±7r/2(m + 
n) , realizing an effective magnetic flux $ = ±7t/2 for \\ / \) g atom and leading to spin-orbit 
coupling and QSH effect 33]. The last term is the Kondo-exchange interaction which 
includes the spin-flip process that scatters the g atom from / i) to |4- / t) by exchanging 
the spin with the \\. / \) e atom at the same site. 
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FIG. 2. (Color online) The cyclotron trajectory of the g atom (a) and the corresponding time evo¬ 
lution of the spin-defect population £ m n m' ri I7m, ra* (t) | 2 (b) at different Kondo-coupling strengths 
V ex = 0 (blue solid), K/4 (red dashed), K (green dotted), and 10 K (purple dash-dotted) in the 
period of Tq. 


III. DAMPED ORBITAL DYNAMICS OF THE MOBILE ATOM 


We revisit the cyclotron dynamics of the single noninteracting | j,) mobile atom in Ref. 30] 
when Kondo-exchange interaction couples the mobile atom to a plaquette of polarized |j") e 
atoms. The motion of the mobile atom here can be derived in the Hilbert space with quantum 

number Stot = Sf ot = § from the state | V(t)) = £ \lm,n(t) c L,n,i + E 


I FM), 
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FIG. 3. (Color online) Time-averaged interaction-induced deviation A £ from the cyclotron orbit 
of noninteracting mobile atom in the period of To as a function of Kondo coupling V ex . And the 
transition from the damped cyclotron orbit to the collapsed trajectory happens at V ex = J + K, 
which is marked by the red star. 


where | FM) denotes the polarized |t) e atoms at four sites A,B,C, and D associated with 
the g atom vacuum 44-49], and S£^ n = n i 7 m,n is the probability amplitude of 

|4) g atom, and that is relevant to the spin-flip process between |4) g atom and |4) e 

atom is the probability amplitude of |4) e spin defect in a background of |4) e atoms. We 
consider the time evolution of the system initially prepared in |T 0 ) = + c^ D \FM) 

with K = 2'kH x (0.27kHz) and J = 2nh x (0.53kHz) in a period of T 0 = 2.3ms as 


Ref. 


30], and during T n , the |4) g atom without Kondo coupling can reproduce the 


loop trajectory in Ref. 


30]. By solving the Schrodinger equation with Eq. (jT]) , we ob¬ 


tain the mean position of the g atom along x and y, ( X) = (Nb + Nc — Na — Nr>)d/2 
and (Y) = (N c + N D — Na — N B )d/2, with the onsite population of the g atom N m ^ n 
((m, n ) G A, R, C, D ) and lattice constant d. The trajectories of the mobile atom for differ¬ 
ent strengths of the Kondo coupling V ex are shown in Fig. [21(a). At weak Kondo coupling 
or around quantum critical point (V ex = K/4,K), the cyclotron orbit gradually shrinks; 
and at strong Kondo coupling (V ex = 1 OK), the loop trajectory is totally destroyed, and 
the atom moves irregularly. Different from cyclotron motion damping induced by repulsive 


interaction in Ref. 


35!] . the damping of the cyclotron orbit of the g atom is attributed 
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to the spin-flip process induced by Kondo-exchange interaction. As a consequence of spin 
flips, the |4,) g atom evolves into a superposition of I'f) and |],), and since the I'f) and |],) 
g atoms manifest opposite chiralities of motion in the spin-dependent magnetic field, the 


cyclotron orbit shall be evidently damped 


30]. Because the spin exchange between |],) g 


atom and |t) e atom creates a |],) spin defect in polarized |t) e atoms, the time evolution 
of the spin-defect population Ylmnm'n' \lm,n i^)\ 2 i 11 Fig. [2(b) monitors the spin-flip event 
during the motion of the g atom. At nontrivial Kondo coupling, the fractional spin defect 
population confirms the occurrence of spin flips and the superposition of I]') and |],) g atoms; 
from V ex = K /4 to K, the spin defect population remarkably increases in the same period 
as a result of strengthening the spin-flip process; and the intriguing oscillation at V ex = 10 K 


can be considered as a signature of the spin-current vortex in Re 


the formation of composite objects under strong spin flips 


49 


, which is related to 


Furthermore, we define 


the time-averaged interaction-induced deviation from the cyclotron orbit of noninteract¬ 
ing |],) mobile atom in the period of T 0 as A £ = jr dt^[AX(t)] 2 + [AY'(f)] 2 , where 
A X(t) = (X(t)) Vex - (X(t)) o and A Y(t) = (Y(t)) Vex - (Y(t))(X(t)) Vex and (Y(t)) Vex 
are the mean position of the g atom along x and y at time t and Kondo coupling strength 
V ex . And a more systematic study of the effect of Kondo coupling on the cyclotron orbit 
is launched by monitoring A£ as a function of V ex in Fig. [3j At weak Kondo coupling 
or around quantum critical point, the deviation Ad increases with the strength of Kondo- 
exchange interaction, indicating that the cyclotron orbit is being further damped, which is 
consistent with the trajectories derived in the same parameter regime in Fig. [2(a). The 
damped cyclotron orbit of the mobile atom begins to collapse when the Kondo coupling V ex 
overwhelms the kinetic energy of the noninteracting mobile atom in the plaquette J + K 
(red star), where the collapsed trajectory manifests itself by the constant deviation Ad with 
varying V ex at strong Kondo coupling. 


IV. COLLECTIVE CYCLOTRON DYNAMICS OF THE KONDO SINGLET 


The damping of the cyclotron orbit above shows that the QSH effect of the g atom 

At strong Kondo coupling, 


30i , [33] is not robust to Kondo-exchange interaction 38 


since the g and e atoms at the same site tend to form a Kondo singlet 2, [49], we naturally 
investigate the behavior of a composite object in the synthetic magnetic field. First we 
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FIG. 4. (Color online) (a) Schematic illustration of the effective magnetic field and diagonal 
tunneling for the Kondo singlet in polarized |t) e atoms, (b) The cyclotron trajectories of the 
Kondo singlet in the period of 2Tq are obtained from H (blue solid), -He// (red dashed) and 
(green dotted). If e spins become nonparallel, for example, by flipping the e spin at B, the cyclotron 
trajectory disappears (purple dash-dotted). 


define the creation operators in a background of |j^) e atoms for the Kondo singlet s' mn — 
^( C L,n,tf!n,n,i ~ and triplets t°^ n = ± (c^f^ + 

and KL = S !^y c L,njL,n,tfrn,n^ where S f m 7 n , with m',n' ± m,n is a restriction by the 
requirement of S tot = Sf ot = | in the Hilbert space. The physical interpretation of this 
definition is that the creation of a composite object is accompanied by the annihilation of 
the e atom at the same site. Then, in the singlet-triplet representation, the Harper-Kondo 











Hamiltonian can be partitioned as 50] 


H = 


V s H t V s 
V t H T V s 


3 Vea 
4 


V a H T V t 


VtH T Vt + 


V e , 


( 2 ) 


V s ,t are the projectors onto the subspaces {|s) = \FM)} and {\t°' 1 ) mn = \FM)} 

whose Kondo-exchange energy are — AM and respectively. Considering that the cou¬ 
pling between {|s) mn } and {0 O,1 ) mn } is almost energetically forbidden by the large Kondo- 
exchange energy gap V ex J, K, the behavior of the Kondo singlet and triplet can be 
considered to be approximately independent. And the orientation of the magnetic held 


can be inferred from ±0, 


of the nearest-neighbor tunneling matrix element — A e ±*^m,n 


m 


{\s) mn } and {0 o,1 ) mri } 33] - In the high-energy triplet subspace, opposite magnetic fields 
simultaneously work during the tunneling of the triplets: \t 0,1 ) mn —> |t°’ 1 ) m+ i n with 
and |t 0,1 ) —> |i 1,0 ) +1 with — 4> m ,n■ Consequently, the triplets’ movement neither real- 


f m,n ' i" /m+l,n 

izes effective magnetic flux nor manifests chirality 30]. However, in the low-energy singlet 
subspace, the Kondo singlet in a background of I]') e atoms experiences a unidirectional 
magnetic held arising from 0 mjri = vr/2(m + n). The effective Hamiltonian in the low-energy 
subspace {|s) mn } is explicitly given by j^l (see Appendix B for more details) 

1 


H eff = V s H t V s + V s H 7 


Ay 

4 V p -‘ 


‘Ptl J v ex 


'y ^ I 


(H T + H K )]P t 

J 


H t V s 


2 e T S m+l,n S m,n + ^ S m,n+l S m,n + h.C. 




_ zj ( s ) I rr( s ) 

- + H d 1 




+ e 


i( f Pn 


)1 ot o I /, /> t 

J ' 5 m+l,n+l' 5 m,n ' n.L.f 


(3) 


where Jd = 4M- and trivial constants are neglected. The hrst term in Eq. @ reconstructs 
the Harper Hamiltonian of Kondo singlet, and the factor 1/2 originates from the fact 
that the effective mass of Kondo singlet is twice the mass of the g atom 49]. The last 
term describes the second-order diagonal tunneling of Kondo singlet from A to C or 
D to B. Then, _|_ e *(<£m,n+</>m, n + i) reflects the interference between two second- 

order hopping paths, which reduces to 1 + e 1 ^ for the paths A —» B —y C and A —)■ 
D —» C. Following the procedure above, the time-reversal counterpart of Eq. (J3J) with 
0m,n —> —0m,n can also be obtained in a background of |]0 e atoms. This means that the 
Kondo singlet experiences a e-spin-dependent magnetic held as a spin-half particle, and the 
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opposite chiralities of the Kondo singlet’s cyclotron orbit shall be locked to polarized ||) 
and ||) e atoms respectively [di], 33], thus realizing the QSH effect of the Kondo singlet as 
a result. From the aforementioned demonstrations where the QSH effect of the g atom is 
destroyed at strong Kondo coupling and the spin-dependent magnetic held on the Kondo 
singlet is irrelevant to the spin state of the g atom, one can realize that the QSH effect of 
the Kondo singlet revealed in our scheme is nontrivial. 

The Harper Hamiltonian and diagonal tunneling for the Kondo singlet in the background 
of ||) e atoms are illustrated in Fig. (U)a). We consider the cyclotron motion of the Kondo 
singlet in the initial state ^j(|s) A + Is)^,) during the period of 2T 0 at V ex = 20/Y In Fig. 
0(b), different trajectories of the Kondo singlet are calculated by H , H e ff and Hy . Clearly, 
the effective Hamiltonian H e ff at strong Kondo coupling resembles the cyclotron dynamics 
obtained by solving the Schrodinger equation with the Harper-Kondo Hamiltonian H and 
accurately describes the cyclotron orbit of the Kondo singlet in the plaquette. Furthermore, 
comparing the cyclotron trajectories obtained by H e ff and H^\ we find an essential role of 
H d in the movement of Kondo singlet. The singlet travels a longer distance under H e ff than 
that predicted by , because H d speeds the singlet up around the plaquette. Moreover, 
the QSH effect of the Kondo singlet is not robust to interactions that destroy the parallel 
alignment of e spins. To verify this statement, we flip one arbitrary spin of the polarized e 
atoms in the initial state, for example, S'^ _ ^(|s) j4 -|-|s) D ), and find that the cyclotron orbital 
motion of the Kondo singlet disappears. The result indicates the competition between the 
QSH effect of the Kondo singlet and the interaction-induced nonparallel state of e spins 42], 


V. EXPERIMENTAL CONSIDERATIONS 


Now we address practical considerations of observing the cyclotron trajectory of the 
Kondo singlet. Here we only consider ultracold 1,1 Yb atoms in the optical lattice. The 
Raman lasers have to be detuned from the resonant frequency of 399 nm of the 1 So— x Pi 
transition. For 1 S 0 atoms trapped in a lattice (A = 532 nm) with the depth V 0 = 15 E R , we 
have tunneling rate t x = 31 Hz and band gap to = 26 kHz. Given Lande factor gj = 5.4 x 1CK 4 
and magnetic field gradient B' ~ 2 x 10 3 mG//im, energy tilting between nearest-neighbor 
sites is estimated to be A ^ 200 Hz. Therefore, the requirement of t x A -C to can be 


satisfied 


51] . The spin-exchange strength V ex on the order of tens of kHz |l7i-il9| can be 
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tuned by adjusting the lattice depth which changes the overlap of local atomic wavefunctions. 
Though the 3 Po atoms can also feel the magnetic field gradient, their wavefunctions are only 
slightly affected because of the large lattice depth. The atom distribution can be detected 


by band-mapping technique in optical superlattices 
be realized by state-of-the-art techniques 


17 


19 


30, 


30]. Key ingredients in our scheme can 


31 


34]. 


VI. CONCLUSION 


In summary, we have studied the cyclotron dynamics of the mobile atom and the Kondo 
singlet in the plaquette with interorbital spin-exchange interaction and predicted a novel 
QSH effect of the spin-zero Kondo singlet in a background of polarized spins, demonstrating 


that the interplay be 
than competition 36- 


ween Kondo-exchange interaction and QSH effect constitutes more 


42] for a composite object. We show that the strong Kondo-exchange 


interaction destroys the loop trajectory of the mobile atom, but instead generates the cy¬ 
clotron orbital motion for a composite object (Kondo singlet). For spin-orbit-coupled heavy 



work may pave a way for experiments on spin-orbit-coupled heavy-fermion systems with 
AEAs in the future. 
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Appendix A: Harper-Kondo Hamiltonian with laser-assisted tunneling 


We extend the sc 
assisted tunneling 


leme of realizing an Abelian gauge field with spin-dependent laser- 


33] to the system considered in the main text as follows: although normal 


tunneling in the x direction is suppressed by a magnetic held gradient A for g atoms with 
opposite magnetic moments, a driven Raman process with driving frequency cu and two- 
photon Rabi frequency hi restores resonant tunneling. We average out rapidly oscillating 
terms in a rotating frame and obtain the time-independent Harper Hamiltonian of the g 
atom: 


Hr = — 




(Ke~ 


ym.n.a. 


C m+l,n,a C m,n,a 


A J C m,n+l,a C m,n,a + h.C .), 


m,n,a 

where = ±0m,n- Subsequently, we generalize the unitary operator in Ref. 

the spin degree of freedom 


(Al) 


51] to 


u = E 


cos(u;t—0 m , n>a )] „t 


C m,n,a C m,n,on 


(A2) 


and derive the expression of Kondo-exchange interaction in the above rotating frame at 
resonant tunneling 


T 

1 f 2 


Hk = V eI V E/t s%„Stf n U + w I ' ti + s‘- tn Si+ n )U 


'-(i> 


„ T 

1 f 2 


iuut ^ 


2 T J t~ dt Ju (7“ Shl ^ m ’ n ) + ^ c ' 


T/ o/,* i 

Kea; / v ° m,n u m,n ' 

m,n 

Vex 'y ] S m,n^m,n A 2 ^"7" S * n ^ S rn,n^m,n A S m,n^m,n)^ (A3) 


where J^(x) are the Bessel functions of the first kind. The anisotropy of the term 
./o( Al sin <5 m: n) can often be relaxed. When ^ < |, because J 0 (^ sin0 mi „) > 0.985, 
the anisotropy for an arbitrary phase 4> m ,n can be reasonably ignored. Thus, in future 
experiments, by choosing a proper ratio between two-photon Rabi frequency and magnetic 
held gradient, the isotropic Kondo-exchange interaction can always be reached in the driven 
optical lattice: 









Appendix B: Projection to derive the effective Hamiltonian (Eq. (3) in the main 
text) 

In the main text, we obtain the effective Hamiltonian in the strong Kondo-coupling regime 
by the projection 

1 


Hrfj — V,II r V s + V s IIr 


V t \-\V ex -(H T + H K )]V t 


HtV s ■ 


Here, we divide H e ff into two parts Hi = V s HtV s and H 2 = H e ff — Hi , and the matrix 
elements of H l 2 are explicitly given in the low-energy singlet subspace {|s) m n } where (m, n) 
represents plaquette sites A, B , C, and D. Hi includes only the nearest-neighbor tunneling 
term 

I S )a I S )b \s)c l S )-D 


f o -K 

(s\Hi\s) mjl = I 


0 


-K o 
2 u 


J 
' 2 


-l \ 

2 

0 


o -{ o -t e ' 


(Bl) 


_ K_ p i7r / 2 

\-j 0 0 j 

which gives Hj, in the main text. H 2 includes the onsite interaction, as well as the nearest- 
neighbor and diagonal tunneling terms. First, the onsite interaction term only contributes 
the trivial constant 


/Inn _ 3 K 2 + J 2 

m,n\S\ n. 2 \S) mn — ^ y 

Second, the nearest-neighbor tunneling term can be safely neglected because 

K 


m+ l,n (s| H 2 |s) m n — m ,n +1 ( s | H 2 |s) m n — 0 


V P , 


~ 0. 


(B2) 


(B3) 


Third and last, we obtain the second-order diagonal tunneling as 

KJ 


m+l,n+l ( s | H 2 \s) m n 


414, 


-(l + e i7r/2 ), 


(B4) 


where 1 + e* 71 '/ 2 originates from the interference between possible paths of diagonal tunneling 
from A to C or D to B. For example, a Hondo singlet initially at site A can pass through 
B or D to arrive at diagonal site C. In the tunneling path A —> B —y C, the singlet carries 
a zero phase because <^4 = 0 and 0^ = 0. However, in the other path A —* D —» C, the 
accumulated phase changes to0 J 4 + 0D = O + | = : |. As a result, the interference between 
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two different tunneling paths contributes e*° + e *2 = 1 + e* 2 . Finally, H -2 gives H^ s> in the 
main text. 
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